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ABSTRACT 


The range uncertainties associated with the refractive atmosphere 
can be mitigated by the technique of two color, or dual wave- 
length, ranging. The precision of the differential time of flight 
(DTOF) measurement depends on the atmospheric dispersion between 
the two wavelengths, the received pulsewidths and photoelectron 
counts, and on the amount of temporal averaging. In general, the 
transmitted wavelengths are not independently chosen but instead 
are generated via nonlinear optics techniques (harmonic crystals, 
Raman scattering, etc.) which also determine their relative pul- 
sewidths. The mean received photoelectrons at each wavelength are 
calculated via the familiar radar link equation which contains 
several wavelength dependent parameters. By collecting the vari- 
ous wavelength dependent terms, one can define a wavelength fig- 
ure of merit for a two color laser ranging system. 

In this paper, we apply the wavelength figure of merit to the 
case of an extremely clear atmosphere and draw several conclu- 
sions regarding the relative merits of fundamental-second har- 
monic, fundamental -third harmonic, second-third harmonic, and 
Raman two color systems. We find that, in spite of the larger 
dispersion between wavelengths, fundamental— third harmonic sys- 
tems have the lowest figure of merit due to a combination of poor 
detector performance at the fundamental and poor atmospheric 
transmission at the third harmonic. Fundamental -second harmonic 
("700 nm and 350 nm) have the highest figure of merit, but 
second-third harmonic systems, using fundamental transmitters 
near 1000 nm, are a close second. Raman— shifted transmitters 
appear to offer no advantage over harmonic systems because of (1) 
the relatively small wavelength separation that can be achieved 
in light gases such as hydrogen and (2) the lack of good ultra- 
short pulse transmitters with an optimum fundamental wavelength 
near 400 nm. 
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1 INTRODUCTION 

With the subcentimeter precisions available from modern satellite 
laser ranging (SLR) hardware [Degnan, 1985], atmospheric refrac- 
tion is a dominant error source in the absolute determination of 
the geometric range from the station to the satellite. While 
atmospheric modelling is believed to reduce the systematic errors 
to roughly one centimeter or less, future progress toward milli- 
meter absolute accuracy ranging will rely on the technique of two 
color, or dual wavelength, ranging. 

In the present paper, we attempt to define optimum wavelengths 
for two color ranging. In order to accomplish this, we must take 
into account all of the wavelength dependent parameters which 
influence our ability to make an accurate differential time of 
flight (DTOF) measurement. As we will see in the ensuing sec- 
tions, a proper accounting of wavelength dependent terms will 
include atmospheric dispersion, atmospheric transmission as a 
function of sea level visibility, transmit antenna and target 
gains, detector responsivities, transmitter availability and pul- 
sewidth, and the detailed characteristics of the available non- 
linear optics techniques for achieving the necessary optical 
frequency translations. 


2 ATMOSPHERIC REFRACTION: THE MARINI-MURRAY MODEL 

In the Marini-Murray model of atmospheric refraction [Marini and 
Murray, 1973], radial variability in the meteorological parame- 
ters (i.e. with altitude) is assumed to be governed by the equa- 
tions for hydrostatic equilibrium, the law of partial pressures, 
and the perfect gas law. This leads to the following equations 
for the spherical range correction, SC^: 
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where \ is the laser wavelength in microns, E is the true eleva- 
tion angle of the satellite in degrees, <f> is the station lati- 
tude, H is the station height above mean sea level, and P H , T H 
and ejj are the surface pressure, temperature, and water vapor 
pressure at the station. The water vapor pressure e H is related 
to the surface percent relative humidity Rh and surface tempera- 
ture Th by the equation 

T H - 273 IS 
237.3 ♦ (T n - 273. 15) 

H ) = YoO 0 ' 1 1 * 1W 

The wavelength dependence of the range correction is contained in 
the dispersion term /(>0 which is plotted in Figure 1. It was 
arbitrarily chosen by Marini and Murray to have a value of unity 
at the ruby laser wavelength of .6943 microns. 
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Figure 1. Atmospheric dispersion in a standard atmosphere as a 
function of wavelength from the near ultraviolet to the near 
infrared. 

3 TWO-COLOR LASER RANGING 

By measuring the pulse times— of— flight at two colors and multi- 
plying the results by the velocity of light in vacuum, c, we 
obtain a measure of the optical path lengths through the 
atmosphere at the two wavelengths. Thus, the atmospheric refrac- 
tion correction is given by 
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(3.1) 


AC = y(L i -L 2 ) = ^j(x l -x 2 ) 


where Li and L 2 are the optical path lengths and and 12 are 
the measured roundtrip times of flight at the two wavelengths 
respectively, and 


rc<7i~ 1 
n g 2 ~ ™g 1 


(3.2) 


where ngi and ng2 are the group refractive indices at the two 
wavelengths. Unfortunately, the wavelength dependence due to the 
"dry" and "wet" components of the atmosphere are different 
[Owens, 1968] whereas the expression for the group refractivity 
Ng used by Marini and Murray assumes no dependence of the water 
vapor term on wavelength, i.e. 

W„(\) = 80.343/(\)£-11.3£ (3.3) 


where P, T, and e are the local pressure, temperature, and water 
vapor partial pressure respectively. Nevertheless, under normal 
conditions of modest humidity, y can be well approximated by the 
expression [Abshire and Gardner, 1985] 

/CM) 

Y /(*> 2 )-/(*-i) (3 ' 4 


where /(>0 is given by (2.2a). If we assume that the two times of 
flight are independently measured, we can express the expected 
variance in the atmospheric correction as 

°^c = (Y c/ 2) 2 (0i + 02) (3.5) 


where and 02 are the RMS errors in the time of flight measure- 
ments at the wavelengths A,j and X 2 respectively., In the ideal 
limit where the differential timing precision is determined only 
by the signal strength, one can write 


0 


AC 


yc 

2 



(3.6) 


where Tpi and r p 2 are the laser pulsewidths and ni and n 2 are the 
received photoelectron signal strengths at the two wavelengths 
respectively. 


4 THE RADAR LINK EQUATION 

The mean signal flux in a range receiver is obtained from the 
familiar radar link equation. The mean number of photoelectrons 
np e recorded by the ranging detector is given by: 
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(4.1) 



where Ti q is the detector quantum efficiency, E-p is the laser 
pulse energy, X is the laser wavelength, h is Planck's constant, 
c is the velocity of light in vacuum, r\ t is the transmit optics 
efficiency, Gt is the transmitter gain, a is the satellite opti- 
cal cross-section, R is the slant range to the target, A r is the 
effective area of the telescope receive aperture, is the 
efficiency of the receive optics, T a is the one-way atmospheric 
transmission, T c is the one way transmissivity of cirrus clouds 
(when present) , and R is the slant range between the station and 
the target. 

In discussing the link equation, we are primarily concerned with 
those terms which exhibit a wavelength dependence. While optical 
coatings certainly exhibit a wavelength dependence thereby 
affecting the transmit and receive optical efficiencies ri t and n r , 
coatings can generally be designed to give approximately equal 
performance once the operating wavelengths are chosen and hence 
will not be included in our discussion. Similarly, experimental 
studies of cirrus cloud transmission have shown no significant 
dependence on wavelength over the band from 0.317 to 12 microns. 
Other terms in (4.1) do have a wavelength dependence which we 
will now discuss. 

4.1 TRANSMITTER GAIN 

A general expression for the transmitter gain is given by 


G«“ 


4n A, 
X 2 




(4.2) 


where A, = na 2 is the area of the transmitting aperture and 

y t ,X) is a geometric factor independent of wavelength 
[Klein and Degnan, 1974]. Note that, for a given transmit 
aperture and a well— collimated system, the transmitter gain is 
inversely proportional to the wavelength squared. 


4.2 TARGET OPTICAL CROSS-SECTION 

The optical cross-section of an unspoiled retroref lector is 
given by [Degnan, 1992] 

«--p (4 ' 3: 

where p is the cube corner reflectivity, A cc = jiR 2 c i s the light 

collecting area of the corner cube, and 4nA cc /X 2 is the on- 
axis retroref lector gain. Even in the presence of complicating 
factors such as velocity aberration and retroreflector 
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spoiling, an array of retroref lectors designed to operate at 
both wavelengths would be expected to retain the same inverse 
square law dependence on wavelength exhibited by (4.3). 

4.3 ATMOSPHERIC ATTENUATION 

In the near-ultraviolet to visible spectral band between 0.3 
and 0.7 n, atmospheric attenuation is dominated by aerosol 
(Mie) scattering but molecular (Rayleigh) and ozone absorption 
also play a role [RCA, 1968]. In the near infrared beyond 
0.7 (x, the plot of atmospheric transmission versus wavelength 
(see Figure 2) is modulated by strong absorption features of 
various molecular constituents in the atmosphere, notably 
water vapor, oxygen, and carbon dioxide. 



Mu) 

Figure 2. Atmospheric transmission as a function of wavelength 
under extremely clear conditions with 2 cm of precipitable 
water vapor at zenith angles of 0, 50, and 70° (corresponding 
to 1, 2, and 3 air masses) respectively. 


The transmission curve presented in Figure 2 corresponds to 
excellent "seeing" conditions (80 Km visibility) and 2 cm of 
precipitable water vapor. It should be noted that atmospheric 
seeing conditions vary widely from day to day and from site to 
site and are usually characterized by "sea level visibility" 
expressed in kilometers. Plots of the sea level attenuation 
coefficient versus wavelength (from 0.4 to 4 n) as a function 
of sea level visibility can be found in the RCA Electro-Optics 
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Handbook 
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TRCA, 19681. For the purpose of this analysis, how 
will consider only the extremely clear atmosphere 
in Figure 2. 


4.4 OPTICAL DETECTORS 

Finally, the availability of high quantum efficiency optical 
detectors at the two laser wavelengths is important. If a com 
mon photocathode is to be used, such as in most s £reak camera 
schemes for performing differential timing, the photocathod 
must be sensitive at both wavelengths. However, since the 
images of the two return pulses can be spatially separated in 
the entrance slit of the streak camera, one can conceive of 
specially constructed streak tubes containing more than one 
photocathode material to obtain the highest sensitivity a 
both wavelengths. It may also be possible, at some future 
date, to do the necessary timing via electronic means, such as 
high speed GaAs technology, without resorting to streak camera 
technology although this capability has not yet been demon 
strated. 


A) “BEST" PHOTOEMISSiVE DETECTOR RESPONSIVITIES 



BESPONSIVITV 


SbhIhhm 

50 

S20 

.35 

60 

K-Cs-Sb - (LIME GLASS) 

.40 

80 

K-Cs-Sb - (LIME GLASS) 

.45 

70 

K-Cs-Sb - (LIME GLASS) 

.50 

63 

K-Cs-Sb + (LIME GLASS) 

.55 

64 

GaAs { - 9741 GLASS) 

.60 

65 

GaAs <- 9741 GLASS) 

65 

67 

GaAs (- 9741 GLASS) 

.70 

68 

GaAs (-*- 9741 GLASS) 

.75 

69 

GaAs ( ^ 9741 GLASS) 

.80 

70 

GaAs ( 9741 GLASS) 

.85 

68 

GaAs ( 9741 GLASS) 

.90 

13 

GaAs ( ■+■ 9741 GLASS) 

.98 

9 

Gain As ( 9741 GLASS) 

1.02 

.8 

GalnAs (+ 9741 GLASS) 

1.1 

.2 

SI 



Figure 3 . Summary of "best" photoemissive detector responsivi- 
ties from the near ultraviolet to the near infrared. 
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Detector sensitivity at a particular wavelength is usually 
expressed as "spectral responsivity" in milliamperes/Watt . It 
is related to quantum efficiency by the equation 

(4. 4 ; 

where is the detector spectral responsivity at wavelength 

K h is Planck's constant, and c is the velocity of light. A 
composite responsivity curve, which is the envelope of indi- 
vidual responsitivity curves for some common visible and near 
infrared photoemissive detectors [Slater, 1980; Zwicker, 19771 
is illustrated in Figure 3. 

5 "OPTIMUM" WAVELENGTHS FOR TWO COLOR SLR 

In choosing "optimum" candidate wavelengths for successful two 
color ranging, there are a variety of technical issues the engi- 
neer must consider. These will be discussed in the ensuing sub- 
sections. As we shall now see, equation (3.6) for the RMS error 
m the atmospheric correction for the photon-limited case, 
combined with the radar link equation (4.1), points the way to 
the selection of a set of optimum wavelengths. Since we want to 
minimize oj^q, the inverse of (3.6) can serve as an overall system 
figure of merit. 

5.1 ATMOSPHERE 

The dependence of (3.6) on the atmospheric dispersive function 
/(A.) illustrates the need for adequate atmospheric dispersion 
between the two wavelengths in order to reduce the severity of 
the timing requirements. The atmospheric dispersion curve in 
1 strongly suggests that one wavelength be chosen to 
lie in the near ultraviolet. On the other hand, atmospheric 
attenuation in the spectral band between 0.3 and 0.7 microns, 
resulting from the combined effects of molecular (Rayleigh) 
and aerosol (Mie) scattering and ozone absorption, also 
increases rapidly in the near ultraviolet as shown in Figure 
2. This will negatively impact the timing precision by lower- 
ing the photoelectron count at the UV wavelength. Furthermore 
in choosing a laser wavelength, it is probably wise to avoid 
the strong water absorption lines in the spectral regions 
between 0.7 and 1.0 microns and beyond 1.1 micron. The high 
variability of water vapor total burden would impact both the 
day-to-day signal strength and cause the pulse group velocity 
to vary via the anomalous dispersion effect near an absorbing 
feature. 

5.2 LASER TRANSMITTER 

The availability of lasers capable of generating high peak 
powers and ultrashort pulsewidths on the order of 35 picosec- 
onds or less is also a consideration. Preference is generally 
given to solid state lasers because of the practical 
difficulties of using liquid dye lasers in the field. Over the 
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past decade, much progress has been made m the developmen 
highly tunable solid state lasers such as Alexandrite (700 to 
810 nm) and Titanium-doped sapphire (600 to 900 nm) . The wid 
bandwidths of these new materials are capable of supporting 
subpicosecond pulsewidths whereas today's workhorse, Nd.YAG, 
is limited to about 10 picoseconds by it]s relatively narrow 
linewidth (120 GHz). However, high bandwidth comes at a price 
- i.e. lower gain - making the construction of high peak power 
Ti: Sapphire and Alexandrite devices more difficult. 


Generally, the wavelengths in two color systems are generated 
from the fundamental wavelength K f via nonlinear optical tech- 
niques such as harmonic generation in crystals or Raman shift- 
ing in gases. This assures simultaneity of emission and 
eliminates (thankfully) the need to synchronize the firings of 
two separate lasers with picosecond precisions. However, 
reliance on nonlinear teqhniques implies that the two wave 
lengths cannot be chosen independently of each other. In the 
case of harmonic generation, the second and third harmonic 
wavelengths are given by 



( 5 . 1 ) 


respectively . 

In Raman shifting, a portion of the incident radiation at 
input frequency, v 0 is shifted by some fixed amount v s (the 
"Stokes shift") toward longer wavelengths relative to the fun- 
damental. One also obtains frequencies at longer ("Stokes") 
and shorter ("Anti-Stokes") wavelengths, but these are 
generally too weak to supply sufficient energy for satellite 
ranging. Large Stokes frequency shifts, and hence high disper 
sion between wavelengths, are obtained by Raman shifting m 
light gases. Hydrogen produces the largest shift of 4155 cm , 
and photon conversion efficiencies as high as 80% have been 
reported. For example, one proposed two color SLR system uses 
the second harmonic of Nd:YAG (532 nm) . in hydrogen to obtain a 
second wavelength output at ^80 nm [Gaignebet et al, 1986]. 

A second consequence of harmonic or Raman generation is that 
the pulsewidth of the secondary wavelength is generally 
shorter than the pulsewidth of the fundamental. From the 
theory of harmonic generation [Degnan, 1979], the harmonic 
pulsewidths are approximately given by 




( 5 . 2 ) 


for low to moderate energy conversion efficiencies (< 50% 
the usual case) . Raman generation depends on third order non— 
linear processes, and the pulsewidth dependence is identical 
to that of third harmonic generation. 
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5.3 WAVELENGTH FIGURE OF MERIT 

In order to treat all potential system configurations on an 
equal basis, some assumptions are in order. We will assume 
that the fundamental laser, from which all other wavelengths 
are derived, is characterized by an energy E and a pulsewidth 
fp which is constant for all wavelengths. Thus, Equation (3.6) 
becomes 


a 


AC 


= ycx p 




(5.3) 


where /3 X and (3 2 are pulsewidth scale factors which depend on 
the nonlinear process used to generate them as in (5.2). In 
addition, we recognize that wavelengths derived via nonlinear 
processes are obtained with some typical energy efficiency 
which we will denote by r|i and r\ 2 respectively. If the funda- 
mental wavelength is used as one of the two wavelengths, we 
will assign values of /9f = 1 and r|/ = 1. For second and third 
harmonic generation in the ultrashort pulse regime, typical 
conversion efficiencies are r) sh = .5 and ti, A =.2 respectively. 

In order to derive a wavelength figure of merit, we must now 
bring together all of the wavelength dependent terms in equa- 
tions (4.1) and (3.5). We obtain for the figure of merit 


F(X 1 ,X 2 ,E)= 


/(M)-/(^ 2 ) 


kt 


Xt 


/(k .) LniP?/?(A. 1 )7'2(\ 1> £') Ti 2 p lR{X 2 }Tl(X 2 ,E) 


(5.4) 

where fix') is the wavelength dispersion term in the Marini- 
Murray atmospheric correction formula, r|, and t ] 2 and and / ? 2 
are the energy conversion efficiencies and pulsewidth 
reduction factors respectively for the relevant nonlinear pro- 
cess, R(X) is the spectral responsivity, and T a {X,E) is the 
one way atmospheric transmission as a function of wavelength 
and elevation angle. The factor of X 4 comes from the combined 
inverse square law dependence of the transmitter and target 
(retroref lector) gains on wavelength. The additional factor of 
X associated with converting detector quantum efficiency to 
spectral responsivity cancels with a similar factor in (4.1) 
which converts transmitter laser energy to the number of 
transmitter photons. In plotting (5.4), we will use the enve- 
lope of the individual photoemitter responsivity curves in 
Figure 3 so that we present each wavelength in its most 
favorable light. 
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6 CONCLUSIONS 

Fiqure 4 provides plots of the wavelength figure of merit as a 
function of the fundamental (laser) wavelength and elevation 
angle for the extremely clear atmosphere depicted in Figure . 
Parts (a), (b) , and (c) correspond to elevation angles of 90, 45, 

and 20 degrees respectively. The three curves within eac P ° 
compare systems which use: (1) the fundamental and second har- 

monic wavelengths; (2) the fundamental and third harmonic wave 
lengths; and (3) the second and third harmonic wavelengths. 

In spite of their greater dispersion, fundamental- third harmonic 
systems have the lowest figure of merit at all elevation angles 
due to a combination of poor detector performance at the funda 
mental and poor atmospheric transmission at the third harmonic. 

The performance of these systems peaks at a fundamental 
wavelength of about .97 microns independent of elevation angle. 

At zenith, fundamental -second harmonic systems, operating at 
wavelengths of 670 and 335 nm, have the highest figure of merit 
(F = 1.75) but second-third harmonic systems, operating at wave 
lengths of 525 and 350 nm (fundamental = 1050 nm) , are almost as 
good (F = 1.6). As one progresses to smaller elevation angles, 
atmospheric attenuation in the ultraviolet begins to dominate and 
the optimum fundamental wavelength is shifted toward longer wave 
lengths with a corresponding reduction in the wavelength figure 
of merit. Thus, the NdrYAG laser, with a fundamental wavelength 
of 1064 nm and a very mature technology, is a near-optimum choice 
for a second-third harmonic system. However, a fundamental second 
harmonic system which utilizes a Ti: Sapphire laser operating in 
the near infrared beyond 670 nm is a possible alternative. 

Figure 5 suggests that a fundamental wavelength of about 400 nm 
is optimum for a hydrogen Raman-shifted laser and that these sys 
terns offer no real advantage over harmonic systems because or 
their lower wavelength figure of merit. This conclusion is 
further supported by the fact that there are no high power solid 
state lasers operating in the near ultraviolet. The principle 
short wavelength devices are excimer ( "excited d imeE ) 9 a s 
charge lasers. At present, excimers cannot achieve ultrashort 
pulsewidths on the order of picoseconds, typically operate in th 
high atmospheric attenuation region of the spectrum below 360 nm, 
and are operationally less desirable than high power solid state 
lasers. Doubling or tripling solid state lasers to achieve a near 
ultraviolet wavelength prior to Raman shifting only decreases the 
overall wavelength figure of merit further by reducing the values 
for the energy efficiency factors rj , and r| 2 . However, this effect 
is partially offset by the slight reduction in pulsewidth (in- 
creased (3 values) resulting from nonlinear generation. 
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(A) E = BO” 



Figure 4. Two color wavelength figure of merit for fundamental- 
second harmonic, fundamental-third harmonic, and second-third 
harmonic systems operating in an extremely clear atmosphere at 
elevation angles of (a) 90°, (b) 45° and (c) 20° respectively. 
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Figure 5. Two color wavelength figure of merit for Raman systems 
operating in an extremely clear atmosphere at elevation angles of 

90°, 45° and 20°. 
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